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ABSTRACT 

Spectropolarimetric observations have been used to map stellar magnetic fields, many of 
which display strong bands of azimuthal fields that are toroidal. A number of explanations 
have been proposed to explain how such fields might be generated though none are defini¬ 
tive. In this paper, we examine the toroidal fields of a sample of 55 stars with magnetic maps, 
with masses in the range 0.1-1.5 M 0 . We find that the energy contained in toroidal fields has 
a power law dependence on the energy contained in poloidal fields. However the power in¬ 
dex is not constant across our sample, with stars less and more massive than 0.5 M G having 
power indices of 0.72+0.08 and 1.25+0.06 respectively. There is some evidence that these 
two power laws correspond to stars in the saturated and unsaturated regimes of the rotation- 
activity relation. Additionally, our sample shows that strong toroidal fields must be generated 
axisymmetrically. The latitudes at which these bands appear depend on the stellar rotation pe¬ 
riod with fast rotators displaying higher latitude bands than slow rotators. The results in this 
paper present new constraints for future dynamo studies. 

Key words: techniques: polarimetric - stars: activity - stars: magnetic field - stars: rotation 


1 INTRODUCTION 


Stellar magnetic fields are ubiquitous across the main sequence, 
playing a crucial role in the evolution of stars. Zeeman-Doppler 
imaging (ZDI), a tomographic imaging technique, has been used 
to map the large-scale surface fields of a wide variety of stars 
such as rapid rotators (|Donati & Collier Cameron||1997| [Donati 


eTaL]fl999 


nati et al. 


|Donati |1999| |Donati et aT~[ 2003| >, M dwarfs jDo 


2008a, Morin et al. |2008b 201 Of , solar-type stars 


(e.g. Petit et ah1 |2008) and planet-hosting stars (e.g. Fares et al. 
|2013 1 . This technique is capable of identifying the field geome¬ 
try as well as the energy budget of the different field components 
(axisymmetric/non-axisymmetric, dipole/quadrupole/higher order 
multipoles, poloidal/toroidal). Additionally, repeated observations 
of individual stars allows the temporal evolution of the field to be 
assessed with magnetic field polarity reversals observed for a num- 


ber of stars (e.g. iDonati, Collier Cameron & Petit 2003; Donati 

et al. 

2008b; Fares et al. 2009; Petit et al. 2009; Morgenthaler et al. 

2011 

i. The wealth of data available has revealed how the mag- 
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netic properties of stars vary with various stellar parameters such 
as Rossby number and age fVidotto et al.|2014) . 

Despite the progress being made in the field of stellar mag¬ 
netism, there are still a number of open questions. For instance, 
what mechanisms are responsible for dynamo action and where 
does the dynamo operate? In the Sun, it is thought that fields are 
mainly generated in a thin layer between the radiative core and 
the convective zone, known as the tachocline, where shearing is 
strongest l jCharbonneau||20K)| . However, the same process can¬ 
not be responsible for field generation in all stars. Below the fully 
convective limit (~0.35M o ), stars do not possess a tachocline and 
hence different dynamo mechanisms must be at work. Addition¬ 
ally, strong toroidal azimuthal fields have been observed on a range 


of stars (e.g. Donati & Collier Cameron 

1997, Donati 1999 Pe- 

tit et al. 2009; Morgenthaler et al. 2012, 

Jeffers et al. 2014, Boro 

Saikia et al. 20151. As highlighted by Jardine et al. ( 2013), such az- 


imuthal bands can only exist if the magnetic field has been stressed 
above its lowest energy state. Additionally, these authors determine 
that stellar winds cannot be the source of these stresses. Other au¬ 
thors have suggested that such azimuthal fields are evidence of dy¬ 
namos distributed t hroughout the convection zone (Donati & CoL 
|lier Cameron|1997]|Donati, Collier Cameron & Petit|2003||Donati 

|et al.|2003|> rather than being confined to the tachocline as in the 
Sun. This view is supported by the work of |Brown et al.| |201 Of 
who demonstrate that a tachocline is not required to generate strong 
bands of toroidal fields in rapidly rotating solar-like stars. 

The toroidal component of stellar magnetic fields contains the 
free energy that, once liberated, is responsible for energetic events. 
Flares, coronal mass ejections and space weather in general have a 


large influence on the stellar environment, and can affect any plan- 

ets orbiting the star (Zarkaj 

2007; Griefimeier, Zarka & Spreeuw 

2007, Llama et a"Li20li]]Vid 

otto et al. 2012, 2013,2015; See et al. 


|2014| [2015. Cohen et al. 20T5). Presently, the toroidal component 
of stellar magnetic fields has only been studied in single, or small 


samples, of stars. Petit et al7[< j2008| ) studied a sample of four solar¬ 
like stars and noted that the stellar rotation period plays an im¬ 
portant role in determining the fraction of magnetic energy in the 
toroidal component of the stellar field. However, the rotation pe¬ 
riod cannot be the sole parameter that determines the toroidal en¬ 
ergy fraction since stars with similar rotation periods show different 
toroidal energy fractions. Additionally, observations of individual 
stars, over multiple epochs, show that the toroidal energy fraction 
can change significantly on the time-s cale of years (DonatijfcClol- 

lier Camero n) 1997 t|Donati et al.| 19991IDonati, Collier Cameron & 

Petitp003jjPetitrtaTJ2009yFarese£alJ201BJJMorgenthalere£al. 

2012||Fares et al.|2013[|Jeffers et al.|2014[|Boro Saikia et al.|2015| >. 

This indicates that the dynamo is dynamic in nature and cannot be 
characterised by single time averaged parameters. 

In this work, we conduct a statistical study of toroidal fields 
using a sample of 55 stars with ZDI maps. We discuss our sample 
of stars and the mathematical representation of their magnetic fields 
in section[2] In section[3]we cover the trends seen within our sample 
with concluding remarks given in section]?] 


2 STELLAR SAMPLE 

For this work, we employ 90 magnetic maps of 55 dwarf stars with 
spectral types spanning F, G, K and M. The original publications 
for each map are referenced in Table [T] Additionally, parameters 
for each star are listed in Table [T] FigTy] shows the rotation peri¬ 
ods and masses of the stars in our sample. Stars with a single map 



Figure 1. The rotation periods and masses of each star in our sample. Filled 
blue points indicate stars observed at one epoch while open red data points 
indicate stars observed over multiple epochs. Dashed lines are included 
showing a stellar mass of 0.5 M 0 and Rossby numbers of 0.01, 0.1 and 
1 . 0 . 


are plotted with filled blue points while stars with multiple maps 
over many epochs are plotted with open red points. It is worth not¬ 
ing that, though we principally focus on mass and rotation rates in 
this work, there may be other parameters, such as stellar age and 
the presence of planets, that affect the magnetic properties of stars. 
Each map contains information about the large-scale magnetic field 
at the surface of the star. As outlined by |Donati et al.|(2006[ >, the 
radial, meridional and azimuthal components of the magnetic field 
of each star (S r ,B 0 ,S 0 respectively) are given by 

B r (6.cf,) = Yj a t."’ Y ^<t>)’ CD 

£,m 

Bg(6, 0 ) = - ^ \Pe,mZe,m(0> 0) + Ve,mXc,m(0, 0 )], (2) 

t,m 

84 ,( 9 , 0) = - ^ {J3e, m Xe, m (9, 0) - ye, m Z{,m(0, 0)], (3) 

e,m 

where 


= c t ,„,PcAcoseW m *, 


(4) 




C(,m dPe Jn (cos6) ^ 

{+1 de e 




C{,m 

£ + 1 


Pf - (C0Sg) »ne^, 
sin 9 


(5) 

(6) 


C{,m — 


2£ + l {£- m)\ 

An (£ + m)\ 


(7) 


In equations 0 - 0 . P[j„{ cos 9) are the associated Legendre poly¬ 
nomials and £ and m are the order and degree respectively. The 
large-scale magnetic field of a star can therefore be fully expressed 
by the a Cm , f} [m and y fm coefficients. These are determined by in¬ 
verting a series of circularly polarised Stokes V profiles [for a more 
detailed discussion on this process, see |Donati et al.| ( ]2006| ]. The 
a (m coefficients characterise the radial component of the field, the 
/3f >m coefficients characterise the meridional and azimuthal com¬ 
ponents of the poloidal field and the y t m coefficients characterise 
the meridional and azimuthal components of the toroidal field. 
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Figure 2. Top: toroidal magnetic energy against poloidal magnetic energy. 
Stars with multiple maps are connected by grey lines. The dotted line in¬ 
dicates (B 2 or ) = (B 2 ol ). The sample is split into stars less massive (pen¬ 
tagon markers) and more massive than 0.5 M G (circle markers). See text 
for further discussion of how these sub-samples were chosen. The two 
dashed lines are best fit lines for these sub-samples; (B 2 or ) °c ($p ol ) a with 
a = 0.72 ± 0.08 and a = 1.25 ± 0.06 for M < 0.5 M© and M > 0.5 M 0 
respectively. AB Dor is shown with triangles. Each point is colour coded by 
stellar mass. Bottom: toroidal energy fraction against stellar mass. Format 
is the same as the top panel. Data-points are also colour coded by stellar 
mass to aid comparison with the top panel. While the two panels show very 
similar information, the difference in behaviours of the two mass ranges is 
much clearer in the bottom panel. 

Throughout the rest of this paper, we will make use of the sur¬ 
face averaged quantity, (B 2 ), where the components of B are given 
by equation 0-0- This is proportional to the average magnetic 
energy density over the surface of a given star. 


3 RESULTS 

3.1 Magnetic energy budget and saturation 

One of the principle ways to split the stellar magnetic energy is into 
its poloidal and toroidal components which is shown in Fig. [2] The 
top panel shows toroidal magnetic energy density against poloidal 
magnetic energy density while the bottom panel shows the toroidal 
energy fraction against stellar mass. Both panels are colour coded 
by stellar mass aiding comparison between them. 

The top panel shows that the toroidal energy is an increasing 
function of poloidal energy. We can attempt to fit a power law of the 
form (B^ or ) oc <Rp 0l >°. However, the sample seems to consist of two 


Figure 3. Poloidal (top) and toroidal (bottom) magnetic energy against 
Rossby number. The formatting is the same as Fig. [2] Right/left facing ar¬ 
rows indicate stars that only have lower/upper estimates for their Rossby 
numbers. The saturated and unsaturated regimes can be clearly seen with 
the transition occurring at a Rossby number of approximately 0.1. Fits 
to the stars in the unsaturated regime, (Bp ol ) Rf 25±019 and (B? or ) °c 
ft- 2 . 99+0.28^ are s | 1()Wn w iUi dashed lines. Note: the magnetic energy axes of 
the two plots are not the same. 


sub-samples. The stars with higher magnetic energy densities ap¬ 
pear to have a smaller power index, a, than the lower energy stars. 
A priori, it is not clear which stars should be included in which 
sub-sample. In the bottom panel, a change of behaviour is evident 
at approximately 0.5 M Q . Stars with a larger mass than this can have 
large toroidal energy fractions but lower mass stars cannot. Though 
the two panels show essentially the same information, this break in 
behaviour at 0.5 M 0 is much clearer in the bottom panel. A num¬ 
ber of authors have previously discussed a sudden change in the 
magnetic properties of M dwarfs at roughly 0.5 M G |Donati et al.j 

|2008a[|Morin et al.|2008bl|20T0||Gregory et al.|2012^ They note 

that this break is roughly coincident with the fully convective limit 
suggesting a link with the change in internal structure. Dividing our 
sample on this basis, we find power index values of a = 0.72 ± 0.08 
and a = 1.25 ± 0.06 for stars less and more massive than 0.5 M G 
respectively. These power laws are plotted in the top panel with 
dashed lines. It is worth noting that, among the M < 0.5 M 0 stars, 
it is the dipole dominate stars that deviate the most from the higher 
index power law. The non-dipolar stars in the bistable regime, as 
discussed by |Morin et al.| ( [2010) l, are roughly compatible with the 
other power law. Additionally, theoretical models predict that these 
non-dipolar stars can vary cyclically and are able to generate sig¬ 
nificant toroidal fields (e.g. Gastine et al.]2013[>. 
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Table 1. Parameters of our sample: star ID, alternative name, stellar mass, rotation period, Rossby number, (B 2 ) (which is proportional to the magnetic energy density), 
toroidal magnetic energy (as a percentage of total energy), axisymmetric magnetic energy (as a percentage of total energy), poloidal axisymmetric magnetic energy 
(as a percentage of poloidal energy), toroidal axisymmetric magnetic energy (as a percentage of toroidal energy) and the observation epoch. Similarly to |Vidotto et al.| 
(2014), we have grouped the stars into solar-like stars, young Suns, hot Jupiter hosts and M dwarfs. References indicate the paper where the magnetic map was originally 
published. For the remaining parameters, references can be found in|Vidotto et al.||2014| and references therein. 


Star 

ID 

Alt. 

name 

M* 

[M G ] 

P rot 

[d] 

Ro 

Solar-like 

stars 

HD 3651 


0.88 

43.4 

1.916 

HD 9986 


1.02 

23 

1.621 

HD 10476 


0.82 

16 

0.576 

HD 20630 

k Ceti 

1.03 

9.3 

0.593 

HD 22049 

e Ei-i 

0.86 

10.3 

0.366 


HD 39587 


HD 56124 

HD 72905 

HD 73350 

HD 75332 

HD 78366 

HD 101501 

HD 131156A 

£ Boo A 

HD 131156B 

f Boo B 

HD 146233 

18 Sco 

HD 166435 

HD 175726 

HD 190771 

HD 201091A 

61 Cyg A 

HD 206860 

HN Peg 


1.03 

4.83 

0.295 

1.03 

18 

1.307 

1 

5 

0.272 

1.04 

12.3 

0.777 

1.21 

4.8 

>1.105 

1.34 

11.4 

>2.781 

0.85 

17.6 

0.663 

0.85 

5.56 

0.256 

0.72 

10.3 

0.611 

0.98 

22.7 

1.324 

1.04 

3.43 

0.259 

1.06 

3.92 

0.272 

0.96 

8.8 

0.453 

0.66 

34.2 

0.786 

1.1 

4.55 

0.388 


Young 

Suns 


Ab Doradus 

0.76 

0.5 

BD-16 351 

0.9 

3.39 

HD 106506 

1.5 

1.39 

HII 296 

0.8 

2.61 

HII 739 

1.08 

2.7 

HIP 12545 

0.58 

4.83 

HIP 76768 

0.61 

3.64 

TYC0486-4943-1 

0.69 

3.75 

TYC5164-567-1 

0.85 

4.71 

TYC6349-0200-1 

0.54 

3.39 

TYC6878-0195-1 

0.65 

5.72 


0.026 


>0.024 


log | j Tor Axi Pol Axi TorAxi Obs 

[% total] [%total] [%pol] [%tor] epoch 


Ref. 


1.19 

3 

87 

87 

98 

- 

-0.29 

7 

53 

50 

94 

- 

0.83 

8 

4 

0 

40 

- 

2.88 

62 

66 

30 

89 

2012 Oct 

2.27 

8 

16 

10 

92 

2007 Jan 

2.08 

6 

58 

59 

45 

2008 Jan 

2.50 

59 

72 

36 

96 

2010 Jan 

2.09 

26 

63 

55 

86 

2011 Oct 

2.63 

45 

59 

43 

80 

2012 Oct 

2.66 

22 

36 

21 

91 

2013 Sep 

2.71 

47 

36 

4 

72 

- 

0.79 

15 

90 

90 

91 

- 

3.03 

82 

81 

13 

97 

- 

2.26 

49 

44 

0 

90 

- 

1.72 

8 

73 

76 

37 

- 

2.35 

4 

91 

93 

29 

- 

2.28 

30 

42 

25 

80 

- 

3.67 

81 

81 

10 

97 

- 

2.60 

32 

42 

24 

81 

- 

0.29 

1 

9 

9 

5 

2007 Aug 

2.73 

36 

49 

37 

70 

- 

2.13 

24 

25 

12 

68 

- 

2.38 

64 

75 

34 

98 

- 

1.05 

4 

59 

58 

96 

- 

2.88 

50 

66 

40 

92 

- 

2.59 

43 

69 

52 

92 

2007 Jul 

2.42 

51 

54 

24 

82 

2008 Aug 

2.21 

11 

45 

48 

28 

2009 Jun 

2.67 

35 

44 

38 

53 

2010 Jul 

2.67 

39 

38 

17 

72 

2011Jul 

2.81 

38 

77 

69 

90 

2013 Jul 


4.92 

28 

18 

9 

41 

2001 Dec 

4.72 

19 

12 

7 

35 

2002 Dec 

3.76 

45 

41 

3 

88 

2012 Sep 

4.51 

53 

47 

17 

75 

2007 Apr 

3.62 

12 

36 

35 

43 

2009 Oct 

2.60 

31 

33 

20 

60 

2009 Oct 

4.35 

35 

60 

47 

83 

2012 Sep 

4.17 

63 

86 

73 

94 

2013 May 

2.99 

23 

25 

12 

67 

2013 Jun 

3.59 

13 

54 

56 

35 

2013 Jun 

3.67 

22 

32 

24 

61 

2013 Jun 

3.68 

30 

38 

19 

83 

2013 Jun 



Petit 

in prep 


Petit 

in prep 

Petit 

in prep 

do Nas 

cime 

nto et al. 

2014 


Jeffers et al. 

5 

4 

Jeffers et al. 

2014 

Jeffers et al. 

2014 

Jeffers et al. 

2014 

Jeffers et al. 

2014 

Jeffers et al. 

20 

4 


Petit 


Petit 


Petit 


Petit 


Petit 


Petit 


Petit 


Petit 


in prep 


l prep 


l prep 


m prep 


l prep 


l prep 


m prep 


in prep 


Petit etaf |20( 8 


Boro Saikia et al 


Boro Saikia et al 


Boro Saikia et al 


Petit 


Petit 


Petit 


Petit 


Petit 


1 prep 


in prep 
in prep 


in prep 


in pre > 


Boro Saikia et al 


Boro Saikia et al 


Boro Saikia et al 


2015 


2015 


(2015 


2015 


2015 


2015 



Donati et al. 

2003 



Donati et al. 

2003 


F 

olsom et al. 

in pre 

p> 


Waite et al. 

20111 

F 

olsom et al. 

in prep 

Folsom et al. 

in prep 

Folsom et al. 

in prep 

Folsom et al. 

in prep 

Folsom et al. 

in prep 

Folsom et al. 

in prep 

Folsom et al. 

in prep 

Folsom et al. 

in prep 


As a brief aside, we discuss the possibility that the smaller 
toroidal energy fraction of stars less massive than 0.5 M 0 is a result 
of the ZDI technique. ZDI captures the large-scale fields but is in¬ 
sensitive to small-scale fields due to flux cancellation effects. [Rein-1 
|ers & ESasri] i j2009[ show that the majority of magnetic flux may be 
missed when the stellar field is reconstructed using only the Stokes 
V signal when compared to the Stokes I signal. If the lowest mass 
M dwarfs have a large fraction of their magnetic energy stored in 


small scale fields, in the form of star spots for example, these fields 
may not be reconstructed by the ZDI technique. However, there is 
no reason to expect ZDI to preferentially miss toroidal field over 
poloidal field. 

It is a well known observational result that the ratio of X-ray to 
bolometric luminosity, R-% = ^x/^bol’ i ncreases with decreasing 
Rossby number, R„, and saturates for stars with Rossby numbers 
less than approximately 0.1 ([Pizzolato et al.||2003[ |Wright et al.| 
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Table 1 - continued 


Star Alt. M* P mt R 0 log | J Tor Axi Pol Axi Tor Axi Obs Ref. 

ID name [M 0 ] [d] [% total] [%total] [%pol] [%tor] epoch 


Hot Jupiter 
hosts 


r Boo 


1.34 

3 

>0.732 

1.00 

63 

63 

17 

90 

2008 Jan 

Fares et al. 

2009 






0.70 

8 

57 

59 

37 

2008 Jun 

Fares et al. 

2009 






0.45 

13 

25 

23 

37 

2008 Jul 

Fares et al. 

2009 






0.85 

12 

33 

35 

15 

2009 May 

Fares et al. 

2013 






1.19 

38 

25 

18 

37 

2010 Jan 

Fares et al. 

2013 






1.02 

31 

38 

28 

59 

2011 Jan 

Fares et al. 

2013 

HD 46375 


0.97 

42 

2.34 

0.70 

1 

77 

77 

86 

2008 Jan 

Fares et al. 

2013 

HD 73256 


1.05 

14 

0.962 

1.74 

22 

20 

3 

79 

2008 Jan 

Fares et al. 

2013 

HD 102195 


0.87 

12.3 

0.473 

2.22 

57 

60 

23 

88 

2008 Jan 

Fares et al. 

2013 

HD 130322 


0.79 

26.1 

0.782 

0.81 

16 

64 

58 

96 

2008 Jan 

Fares et al. 

2013 

HD 179949 


1.21 

7.6 

>1.726 

0.83 

18 

58 

54 

78 

2007 Jun 

Fares et al. 

2012 






1.13 

8 

34 

34 

32 

2009 Sep 

Fares et al. 

2012 

HD 189733 


0.82 

12.5 

0.403 

2.69 

57 

62 

26 

91 

2007 Jun 

Fares et al. 

2010 

M dwarf 

stars 





3.10 

76 

77 

16 

96 

2008 Jul 

Fares et al. 

2010 

GJ 569A 

CE Boo 

0.48 

14.7 

<0.288 

4.17 

5 

96 

96 

100 

2008 Jan 

Donati et al. 

2008a 

GJ 410 

DS Leo 

0.58 

14 

<0.267 

4.10 

82 

92 

58 

99 

2007 Jan 

Donati et al. 

2008a 






4.02 

80 

78 

15 

94 

2007 Dec 

Donati et al. 

2008a 

GJ 182 


0.75 

4.35 

0.054 

4.60 

68 

66 

15 

90 

2007 Jan 

Donati et al. 

2008a 

GJ 49 


0.57 

18.6 

<0.352 

2.94 

52 

84 

67 

100 

2007 Jul 

Donati et al. 

2008a 

GJ 388 

AD Leo 

0.42 

2.24 

0.047 

4.77 

1 

96 

97 

22 

2007 Jan 

Morin et al. 

2008b 






4.75 

3 

89 

92 

8 

2008 Jan 

Morin et al. 

2008b 

GJ 494A 

DT Vir 

0.59 

2.85 

0.092 

4.44 

62 

61 

12 

91 

2007 Jan 

Donati et al. 

2008a 






4.51 

47 

49 

17 

85 

2007 Dec 

Donati et al. 

2008a 

GJ 896 A 

EQ Peg A 

0.39 

1.06 

0.02 

5.31 

12 

66 

71 

29 

2006 Aug 

Morin et al. 

2008b 

GJ 896 B 

EQ Peg B 

0.25 

0.4 

0.005 

5.35 

2 

92 

93 

42 

2006 Aug 

Morin et al. 

2008b 

GJ 873 

EVLac 

0.32 

4.37 

0.068 

5.61 

8 

31 

28 

60 

2006 Aug 

Morin et al. 

2008b 






5.48 

2 

29 

29 

18 

2007 July 

Morin et al. 

2008b 

GJ 1111 

DXCnc 

0.1 

0.46 

0.005 

4.14 

5 

77 

77 

65 

2007 

Morin et al. 

2010 






3.85 

19 

40 

30 

82 

2008 

Morin et al. 

2010 






3.83 

28 

65 

61 

77 

2009 

Morin et al. 

2010 

GJ 1156 


0.14 

0.49 

0.005 

3.71 

8 

3 

2 

18 

2007 

Morin et al. 

2010 






4.21 

11 

17 

12 

57 

2008 

Morin et al. 

2010 






4.05 

5 

1 

1 

7 

2009 

Morin et al. 

2010 

GJ 1245B 


0.12 

0.71 

0.007 

4.59 

14 

10 

6 

37 

2006 

Morin et al. 

2010 






3.72 

9 

15 

13 

30 

2008 

Morin et al. 

2010 

GJ 9520 

OTSer 

0.55 

3.4 

0.097 

4.38 

33 

70 

61 

90 

2008 Feb 

Donati et al. 

2008a 

V374 Peg 


0.28 

0.45 

0.006 

5.75 

3 

79 

81 

5 

2005 Aug 

Morin et al. 

2008a 






5.61 

3 

76 

79 

2 

2006 Aug 

Morin et al. 

2008a 

GJ412B 

WXUMa 

0.1 

0.78 

0.008 

6.16 

2 

92 

92 

89 

2006 

Morin et al. 

2010 






6.35 

2 

92 

94 

38 

2007 

Morin et al. 

2010 






6.34 

2 

84 

86 

3 

2008 

Morin et al. 

2010 






6.56 

3 

94 

96 

20 

2009 

Morin et al. 

2010 

GJ 285 

YZCmi 

0.32 

2.77 

0.04 

5.72 

5 

57 

59 

19 

2007 Jan 

Morin et al. 

2008b 






5.66 

2 

87 

88 

25 

2008 Jan 

Morin et al. 

2008b 


|2011| . Here, the Rossby number is defined as the ratio of the stel¬ 
lar rotation period to the convective turnover time. |Vidotto et al.| 
l |2014j l showed that stellar magnetism has the same qualitative de¬ 
pendence on Rossby number. However, this result was derived us¬ 
ing the radial component of the surface fields only, and hence does 
not consider the toroidal field. In Fig. [3] we plot the poloidal (top 
panel) and toroidal (bottom panel) magnetic energy densities as a 
function of Rossby number. We use the same Rossby number es¬ 
timates as [Vidotto et al.| < [2014] >, where further discussion of the 
estimates can be found. Left and right facing arrows indicate stars 


where the Rossby number estimate is only an upper or lower limit 
respectively. 

It is clear to see that both components qualitatively follow the 
well known behaviour with some quantitative differences. In both 
cases, the cutoff between the saturated and unsaturated regimes oc¬ 
curs at R„ ~ 0.1. However, in the saturated regime, the average 
magnetic energy of the poloidal fields is higher than that of the 
toroidal fields by just over an order of magnitude. Additionally, in 
the unsaturated regime, the slope is steeper for the toroidal com¬ 
ponent. The similar behaviour indicates that the same mechanism 
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is responsible for generating both components or that both compo¬ 
nents are generated from each other. 

|Wright et al.| |20TTJ suggested that the difference in R x be¬ 
haviour in the saturated and unsaturated regimes can be attributed 
to different dynamo mechanisms operating in each of the regimes 
rather than any actual saturation effect. They argued this on the 
basis that the age at which stars transition from the rotational C se¬ 
quence to the I sequence l jBarnes|2(X)3ft is coincident with the tran¬ 
sition from saturated to unsaturated regimes at R„ ~ 0.1 (see their 
Fig. 4). It is interesting to note that the majority of the stars in the 
saturated regime in Fig. [3] have M < 0.5 M 0 (pentagon symbols). 
Assuming that the unsaturated stars are those more massive than 
0.5 M 0 , we find fits of (B 2 ^) oc R ~ 2 - 25±019 and (B 2 m ) oc R ~ 2 - 99±0 ' 28 
for the poloidal and toroidal fields respectively. These fits are plot¬ 
ted with dashed lines. It is worth noting that these two fits and the fit 
in Fig.0with the higher power index, a = 1.25, are the three possi¬ 
ble 2D projections of the same relatively tight sequence of stars in 
, R 0 ) parameter space. 

If there is a physical basis for the saturated stars having masses 
less than 0.5 M G , the fact that they appear to have a different power 
index, a, to the rest of the stars, in Fig. [2] is further evidence for the 
suggestion of |Wright et al.[ i 2011| l. However, care must be taken 
with this interpretation. Fig. I] shows that the majority of stars in 
our sample that are less than 0.5 M 0 also have Rossby numbers 
less than 0.1. Therefore, our interpretation is not that stars in the 
saturated regime must be less massive than 0.5 M G . Rather, it is 
that stars in the saturated regime are those with Rossby numbers 
less than ~0.1 which also happen to be stars less massive then 
0.5 M 0 in our sample. Noticeably, with the exception of AB Dor 
(R„ = 0.028, M = l.OM 0 ), there is a dearth of stars with Rossby 
numbers less than ~0.1 and masses bigger than 0.5 M 0 in our sam¬ 
ple. That the lowest mass stars in our sample are also the fastest 
rotators is not surprising since these stars spin down less rapidly 
than higher mass stars [Mohanty & Basri|2003[ . Correspondingly 
R„ < 0.1, M* > 0.5 M 0 stars are harder to find and map. Indeed a 
number of authors have previously commented on the difficulty in 
separating stellar mass and Rossby number effects due to this bias 
( |Donati et al.||2008a| |Morin et al.||2008bl |2010[ |Reiners & Basri| 


|2009| Gastine et al.|2013| >. Even when we consider AB Dor (plotted 
with triangles in Figs. [2] and [3}, which does fall into the R„ < 0.1, 
M, > 0.5 M 0 region of parameter space, the picture does not be¬ 
come any clearer. In the top panel of Fig. [2] AB Dor lies relatively 
close to both sub-samples such that it is difficult to tell which sub¬ 
sample it would be more appropriate for it to be in. Likewise, in 
Fig. 0 AB Dor falls in an intermediate region. On one hand, one 
might consider it a saturated star on the basis of its Rossby number 
which is smaller than R a ~ 0.1, the value typically used to de¬ 
lineate the saturated and unsaturated regimes. However, one might 
also consider it to be an unsaturated star on the basis that it follows 
the same trend as the other unsaturated stars, lying at the tail end of 
that sequence. If the hypothesis of [Wright et al.|(201 1 [ is correct, 
one might expect the magnetic fields of the stars in this region of 
parameter space to obey the power law with the smaller index from 
Fig-0 

A potential problem with this interpretation is that our two 
sub-samples were initially chosen on the basis of stellar mass, 
since this parameter better discriminates between the two power 
laws in Fig. 0 rather than Rossby number. It may be the case that 
the two power laws in Fig. 0 do not correspond to the saturated 
and unsaturated regimes. Indeed, until more stars in the R a <0.1, 
M* > 0.5 M e region of parameter space (corresponding to stars 
of rotation periods of less than a few days - see Fig. 01 have their 
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Figure 4. The latitude at which the toroidal azimuthal field peaks for each 
star as a function of stellar rotation period. Each point is colour coded by 
stellar mass. The size of the points indicates the full-width half-maximum 
value for the band (see inset). Stars with multiple maps are connected by 
grey lines. Numerical values for the data points are included in Table0 


surface fields mapped, it will be difficult to conclusively confirm or 
reject this hypothesis. 

Recently, Reiners, Schiissler & Passegger|[2014) reinterpreted 
the data of |Wrightet al.[[201 I [ These authors show that, in the un¬ 
saturated regime, Rx = Px/^bol shows less scatter when plotted 
against R~*P~ 2 t , where R, and P rot are the stellar radius and rotation 
period respectively, rather than Rossby number, as is traditional in 
these types of studies. This formulation is approximately equiva¬ 
lent to L-^ K Pmv i e. that X-ray luminosity depends only on the 
stellar rotation period, in the unsaturated regime. The authors make 
no claims as to whether Rx as a function of R a or Lx as a function 
of P rot is the more physically fundamental relationship. When plot¬ 
ting (Rp 0l ) and (B 2 m ) against R~*P~ 2 t (not shown), we once again 
find that the data separates into saturated and unsaturated regimes. 
Similarly to Fig. 0 the stars in the saturated regime are those that 
are less massive than 0.5 M 0 . The suggestion that different dynamo 
mechanisms are present in the saturated and unsaturated regimes is 
therefore possible under either interpretation of the rotation-activity 
relation. 


3.2 Band latitudes 

In this section we determine the latitude and extent of azimuthal 
bands of toroidal field in our sample of stars. Firstly, by examin¬ 
ing the magnetic maps, we eliminate those stars that show no evi¬ 
dence of strong bands by determining the fraction of latitudes with 
only a single polarity in the toroidal component of the azimuthal 
field. From a visual inspection, stars in which this fraction is less 
than 0.2 do not have clear bands and we do not attempt to find a 
band latitude for these stars. For the remaining 67 maps, we av¬ 
erage the field strengths over every longitude to obtain an average 
field strength as a function of latitude and then take the absolute 
value, i.e. \{B tor ^)\(X) where A is the latitude. We plot the latitude 
at which this function is maximal in Fig.0as a function of rotation 
period for each star. Stars with maps over multiple epochs are in¬ 
dicated by grey lines. An indication of the band width is given by 
the full-width half-maximum of the peak in \{B, or( f,)\{A). The size of 
each data point is scaled by the full-width half maximum. All points 
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Table 2. Numerical values for Fig.^] For each star the rotation period, latitude of the band and the full-width half-maximum of the band are 
listed. 


Star 

ID 

P rot 

[d] 

Latitude 

[degrees] 

FWHM 

[degrees] 

Star 

ID 

P rot 

[d] 

Latitude 

[degrees] 

FWHM 

[degrees] 

HD 3651 

43.4 

16.9 

47.8 

TYC5164-567-1 

4.71 

64.7 

21.1 

HD 9986 

23 

22.5 

47.8 

TYC6349-0200-1 

3.39 

73.1 

14.1 

HD 10476 

16 

47.8 

35.2 

r Boo (2009 May) 

3 

64.7 

25.3 

HD 20630 

9.3 

25.3 

40.8 

r boo (2010 Jan) 

3 

30.9 

25.3 

HD 22049 12007 Jan) 

10.3 

42.2 

42.2 

HD 46375 

42 

33.8 

42.2 

HD 22049 (2010 Jan) 

10.3 

33.8 

49.2 

HD 73256 

14 

16.9 

47.8 

HD 22049 (2011 Oct) 

10.3 

36.6 

42.2 

HD 102195 

12.3 

39.4 

42.2 

HD 22049 (2012 Oct) 

10.3 

33.8 

49.2 

HD 130322 

26.1 

16.9 

49.2 

HD 22049 (2013 Oct) 

10.3 

36.6 

46.4 

HD 179949 (2009 Sep) 

7.6 

59.1 

25.3 

HD 39587 

4.83 

22.5 

42.2 

HD 189733 (2007 Jun) 

12.5 

2.8 

43.6 

HD 56124 

18 

47.8 

35.2 

HD 189733 (2008 Jul) 

12.5 

5.6 

52 

HD 72905 

5 

25.3 

40.8 

GJ 569A 

14.7 

30.9 

39.4 

HD 73350 

12.3 

14.1 

45 

GJ 410 (2007 Jan) 

14 

28.1 

47.8 

HD 78366 

11.4 

53.4 

42.2 

GJ 410 (2007 Dec) 

14 

28.1 

46.4 

HD 101501 

17.6 

16.9 

47.8 

GJ 182 

4.35 

14.1 

32.3 

HD 131156A 

5.56 

36.6 

47.8 

GJ 49 

18.6 

45 

38 

HD 131156B 

10.3 

47.8 

35.2 

GJ 388(2007 Jan) 

2.24 

36.6 

32.3 

HD 166435 

3.43 

16.9 

32.3 

GJ 494A (2007 Jan) 

2.85 

25.3 

45 

HD 175726 

3.92 

59.1 

22.5 

GJ 494A (2007 Dec) 

2.85 

42.2 

47.8 

HD 190771 

8.8 

33.8 

39.4 

GJ 896 A 

1.06 

59.1 

23.9 

HD 201091 

34.2 

30.9 

45 

GJ 896 B 

0.4 

75.9 

12.7 

HD 206860 

4.55 

25.3 

39.4 

GJ 873 (2007 July) 

4.37 

53.4 

15.5 

HD 206860 (2007 Jul) 

4.6 

36.6 

57.7 

GJ 1111 (2007) 

0.46 

19.7 

28.1 

HD 206860 (2008 Aug) 

4.6 

8.4 

32.3 

GJ 1111 (2008) 

0.46 

19.7 

42.2 

HD 206860 (2010 Jul) 

4.6 

8.4 

25.3 

GJ 1156(2007) 

0.49 

25.3 

30.9 

HD 206860 (2011 Jul) 

4.6 

8.4 

53.4 

GJ 1156 (2009) 

0.49 

39.4 

16.9 

HD 206860 (2013 Jul) 

4.6 

14.1 

50.6 

GJ 9520 

3.4 

42.2 

30.9 

HD 106506 

1.39 

61.9 

19.7 

GJ 1245B (2006) 

0.71 

14.1 

25.3 

HII 296 

2.61 

16.9 

23.9 

GJ 1245B(2008) 

0.71 

45 

38 

HII 739 

2.7 

56.2 

18.3 

GJ412B (2006) 

0.78 

25.3 

28.1 

HIP 12545 

4.83 

59.1 

38 

GJ412B (2009) 

0.78 

19.7 

19.7 

HIP 76768 

3.64 

64.7 

25.3 

GJ 285 (2008 Jan) 

2.77 

33.8 

38 

TYC0486-4943-1 

3.75 

5.6 

42.2 






are colour coded by stellar mass. Numerical values for this plot are 
included in Table[2] It is worth highlighting that multiple azimuthal 
bands of opposing polarity were seen in the rapid rotators HR 1099 
and AB Dor l |Donati, Collier Cameron & Petit|2003| . 

Several authors have noted that magnetic flux tends to emerge 
at higher latitudes on stars with shorter rotation periods as a result 
of a larger Coriolis force dominating over the magnetic buoyancy of 
the flux tubes ( [Schuessler & Solanki|1992{ |Schuessler et al.|1996| 

|Granzer|2002[ l. Our results are in qualitative agreement of this state¬ 

ment, with the upper envelope of band latitudes in Fig.[4]showing a 
decreasing trend at longer rotation periods. However, the interplay 
between Coriolis and buoyancy forces alone cannot explain the 
large range of band latitudes seen at a given rotation period. This 
is especially true of stars with multiple maps that show the band 
latitude changing significantly over the course of months/years. In 
some cases, the azimuthal bands even disappear and reappear be¬ 
tween observation epochs. 


3.3 Field orientation 

From a visual inspection of the magnetic maps, toroidal azimuthal 
bands appear predominantly axisymmetrically (e.g. |Jardine et al.| 
|2013[ >. Figure [5] shows the fraction of magnetic energy in axisym- 
metric modes against the fraction contained in the toroidal compo¬ 


nent of the field. Mathematically, these correspond to the energy in 
the m = 0 modeQand the y^„, coefficients respectively. Stars with 
multiple epochs of observations are connected by grey lines and 
each point is colour coded by stellar mass. The most striking fea¬ 
ture is the dearth of toroidal non-axisymmetric stars, in the lower 
right corner of the plot. Additionally, there is no trend evident with 
stellar mass. 

The striking trend in Fig. [5] warrants a check for potential bi¬ 
ases in the ZDI technique. Crosstalk between the field components, 
especially the radial and meridional fields is well known and has 
been characterised iDonati & Brown 1997^ In particular, we are 
interested in energy from the toroidal field leaking into the poloidal 
field. To this end, we conducted a systematic series of tests. Specif¬ 
ically, we created a grid of synthetic Stokes profiles corresponding 
to magnetic field geometries where the axisymmetric and toroidal 


1 It is worth noting that an alternative definition of axisymmetric fields, 
m < if/2, is sometimes used in the literature. This definition encapsulates 
modes that are nearly aligned with the rotation axis whereas the m = 0 
definition only accounts for modes that are exa ctly aligned. 

2 In the tests presented by Donati & Brown l | 1 997f , the individual field 
components are considered to be independent parameters. With the new 
spherical harmonic implementation, as described by [Donati et al.|{2006} , 
crosstalk is considerably reduced. 
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Figure 5. Percentage of total magnetic energy contained in the axisym- 
metric component of the field against percentage of total magnetic energy 
contained in toroidal field. Stars observed at multiple epochs are connected 
by grey lines. The data points are colour coded by stellar mass and scale 
with the total magnetic energy density, (B 1 ) (see inset). [Rosen , Kochukhov] 
|& Wade| ^2015) reconstructed the surface field of II Peg using Stokes IV 
(circularly polarised light) and also Stokes IQUV (circularly and linearly 
polarised light). These data points are shown by black diamond and square 
markers respectively. 


energy fractions ranged from zero to one. Magnetic maps were 
then reconstructed from these profiles and compared to the origi¬ 
nal maps that the synthetic profiles were created from. Full details 
of these tests can be found in Folsom et al. (in prep b). We find that 
there is some energy leakage from toroidal into poloidal fields (and 
vice versa). This effect is largest for geometries where the order, f, 
is equal to the degree, m and smallest when there is a large differ¬ 
ence between ( and m. Flowever, the effect is not large enough to 
explain the dearth of points in the lower right hand corner of Fig. [5] 

The magnetic maps used in this study were all reconstructed 
using the Stokes I (unpolarised) and Stokes V (circularly polarised) 
profiles with one exception. |Rosen, Kochukhov & Wade] ( |2015) ob¬ 
served the star II Peg, during two epochs, in all four Stokes pa¬ 
rameters (IQUV; unpolarised, linearly polarised and circularly po¬ 
larised). They subsequently reconstructed magnetic maps using the 
more commonly used two Stokes profiles (IV) and again using all 
four. II Peg is plotted on Fig.[5]with black squares and diamonds in¬ 
dicating results obtained using Stokes IQUV and Stokes IV respec¬ 
tively. Rosen, Kochukhov & Wade ^2015) found that the toroidal 
and axisymmetric energy fractions are similar for both reconstruc¬ 
tions though the maps constructed with Stokes IQUV contained 
significantly more energy than the maps constructed from Stokes 
IV. The latter effect is unsurprising since one should expect more 
information to be reconstructed, and hence more fields, when using 
more data. The data points of Rosen, Kochukhov & Wade (2015| 
fall within the trend shown by the rest of the sample. II Peg follows 
this trend even though it is evolving off the main sequence whereas 
the rest of the sample are less evolved. 

Further insights may be gained by splitting the axisymmetric 
energy into its poloidal and toroidal components. The top panel of 
Fig. 0 shows the fraction of poloidal energy that is axisymmetric, 
i.e. E(ae, m =o,l3e,m=o)/E(a e , m ,/3e,m), against the toroidal energy frac¬ 
tion. The bottom panel of Fig. [6] is similar but plots the fraction 
of toroidal energy that is axisymmetric, i.e. E(y liln=0 )l E(y ej ^). The 
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Figure 6. Axisymmetric poloidal energy as a percentage of total poloidal 
energy (top panel) and axisymmetric toroidal energy as a percentage of total 
toroidal energy (bottom panel) as a function of toroidal energy. Format is 
the same as Fig. [5] 


formatting is the same as that of Fig. 5] As before, the data points 
of |Rosen, Kochukhov & Wade| ( [2015 lfbr II Peg are in agreement 
with the rest of the sample. It is clear that the poloidal and toroidal 
fields behave differently. While the axisymmetric poloidal energy 
does not show a clear trend, a clear one is present for the axisym¬ 
metric toroidal energy. The data show that toroidal fields are gener¬ 
ated in a preferentially axisymmetric manner. This suggests that the 
toroidal field generation mechanism is sensitive to the rotation axis 
in a way that the poloidal field is not. Noticeably, there is a cluster 
of M dwarfs which have dipole dominated fields in the top left hand 
corner of the plot (these are also present in the top left of Fig. [5}. 
|Gastine et al.|(2013) proposed that the strong dipolar component of 
these stars inhibits differential rotation and hence the generation of 
strong toroidal fields through the Omega effect. 

We can also consider how the poloidal and toroidal fields are 
oriented with respect to each other. Figure [7] plots the fraction of 
poloidal energy that is axisymmetric against the fraction of toroidal 
energy that is axisymmetric, i.e. the ordinates of both panels in Fig. 
[6] The data points fill the entire available parameter space for this 
plot indicating that the poloidal and toroidal field orientations are 
not constrained by each other. Alternatively, their orientations may 
be related to each other but in a complex and non trivial way. It 
is particularly interesting that while the orientations of the poloidal 
and toroidal fields are uncorrelated, there is a relatively tight depen¬ 
dence of the toroidal energy density on the poloidal energy density 
(see Fig. [2]. 
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Figure 7. Axisymmetric poloidal energy as a percentage of total poloidal 
energy against axisymmetric toroidal energy as a percentage of total 
toroidal energy. Format is the same as Fig. [5] Similarly to Fig. [5] the in¬ 
set shows (B 2 ) values of 5 X 10 _1 G 2 to 5 X 10“G 2 in multiples of 10 for the 
smallest to biggest points. The numerical values are omitted due to lack of 
space. 


4 CONCLUSION 

Over the last decade, large strides have been made in characteris¬ 
ing the strength, geometry and time evolution of cool star magnetic 
fields. Many stars display strong toroidal azimuthal bands, though 
currently there is no definitive explanation for their formation. In 
this paper, we explore the strength and geometry of toroidal mag¬ 
netic fields by analysing 55 stars that have had their surface mag¬ 
netic fields mapped with ZDI. This is the first time toroidal fields 
have been studied in a large sample of stars. The results presented 
show clear constraints for future dynamo models and can be sum¬ 
marised as follows. 

Our sample shows that strong toroidal fields are strongly ax¬ 
isymmetric. Conversely, the orientation of poloidal fields are un¬ 
constrained with respect to the rotation axis. There is also evidence 
that the underlying dynamo type affects the relative strength of the 
poloidal and toroidal fields. We find that stars less massive than 
0.5 M 0 tend to have small toroidal energy fractions while more 
massive stars can have substantially higher toroidal energy frac¬ 
tions. Additionally, the toroidal field energies of the M > 0.5 M 0 
stars have a steeper power law dependence on the poloidal field en¬ 
ergy than those of the M < 0.5 M 0 stars. These results are in line 
with the results of |Donati et al.| l |2008a| >; |Morin et al.| ( |2008b[[2010[ >; 
Gregory et al.] p012| who noted changes in magnetic properties 
at 0.5 M g , roughly coinciding with the fully convective limit. This 
suggests some link between the break in magnetic properties with 
the change in internal structure. These authors also note a strong 
break in the amount of differential rotation at 0.5 M 0 . It is perhaps 
telling that the stars with the strongest toroidal fields are those that 
display differential rotation at their surface. On the other hand, the 
lowest mass stars, which do not have large toroidal energy fractions 
(see Fig. [2]), display little to no differential rotation. This may be a 
hint that an Omega effect, in the form of differential rotation, is re¬ 
sponsible for generating axisymmetric toroidal fields. It would be 
interesting to compare these observational results to the theoreti¬ 
cal dynamo simulations of|Brown et al.| ([20 10[t and see if their ax- 


isymmetry and poloidal/toroidal energy fractions follow the same 
trends. 

We find that the latitude of the azimuthal bands depends, in 
part, on the stellar rotation period. On fast rotators, buoyant flux 
tubes feel a strong Coriolis force and are deflected polewards al¬ 
lowing flux emergence to occur at higher latitudes. Additionally, 
we find bands at low to mid latitudes for all rotation periods. There¬ 
fore, rotation period cannot be the only parameter determining the 
latitude of flux emergence. However, at present, it is unclear what 
further parameters may be important. 

Additionally, we find that the poloidal and toroidal energies 
both follow the same qualitative behaviour with Rossby number. 
Both show saturation at low Rossby number with magnetic ener¬ 
gies decreasing at higher Rossby number. This is also the qualita¬ 
tive behaviour that X-ray luminosity normalised to bolometric lu¬ 
minosity shows as a function of Rossby number, baring some quan¬ 
titative differences such as the slope in the unsaturated regime. That 
the two field components behave similarly indicates that they share 
a common generation mechanism or that they are generated from 
each other. There is evidence that, in our sample, the M < 0.5 M 0 
stars correspond to stars in the saturated regime giving credence to 
the claim of |Wright et al.|f20 1 I ^ that stars in the saturated and un¬ 
saturated regimes posses different dynamo mechanisms. However, 
it remains to be seen whether this is an effect of the bias in our 
sample. Future spectropolarimetric observations of cool stars with 
R a <0.1, M* > 0.5 M g will help confirm or disprove this hypothe¬ 
sis. 
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